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Abstract 

We investigate a hybrid structure consisting of 20±4 implanted 31 P atoms close to a gate-induced 
silicon single electron transistor (SiSET). In this configuration, the SiSET is extremely sensitive 
to the charge state of the nearby centers, turning from the off state to the conducting state when 
the charge configuration is changed. We present a method to measure fast electron tunnel rates 
between donors and the SiSET island, using a pulsed voltage scheme and low-bandwidth current 
detection. The experimental findings are quantitatively discussed using a rate equation model, 
enabling the extraction of the capture and emission rates. 
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The readout of a single spin is one of the key elements in spin-based quantum information 
processing schemes [HI2]- One may distinguish between single-shot readout, where the pro- 
jective measurement of a single spin is performed in real time, and "spectroscopic" readout, 
where the expectation value of the spin state is deduced from a time-averaged quantity (e.g. 
electrical current, fluorescence emission, ...). Single-shot readout has been demonstrated 
in GaAs/AlGaAs quantum dots j3j 0], while spectroscopic readout has been obtained in a 
variety of systems, from quantum dots [5], [8j [9] to NV centers in diamond [10J and 
dopant atoms in silicon [TT]. Carbon and silicon are particularly attractive platforms for 
solid-state spin-based quantum processors, because they can be isotopically purified to min- 
imize decoherence induced by nuclear spins. However, single-shot readout in these systems 
has not yet been demonstrated. Recently, an architecture for single-shot readout of a donor 
spin in silicon was proposed [12], consisting of a single implanted P donor [13] in close prox- 
imity to an induced silicon single-electron transistor (SiSET) [13] . The approach employs 
a readout principle similar to the one successfully demonstrated in GaAs/AlGaAs single 
quantum dots [3], where the spin state of the electron is deduced from the time-resolved 
observation of spin-dependent tunneling between the dot and a charge reservoir. However, in 
the donor-based proposal [12], the bulk charge reservoir is replaced by the island of a SiSET. 
This configuration is predicted to yield very large charge transfer signals, thereby allowing 
high-fidelity single-shot spin readout. The timescale of the projective spin measurement is 
set by the electron tunneling time between donor and SiSET, which must be controlled and 
understood before attempting spin readout. 

In this Letter we demonstrate and investigate the tunneling of electrons in a hybrid device, 
consisting of approx. 20 ±4 31 P donors, implanted next to an induced SiSET. We show that 
the current through the SiSET, /set, can be switched from zero to the maximum value by 
transferring an electron from a charge center to the SET island. By applying voltage pulses 
to a gate near the donors while monitoring /set, we study the probability for an electron to 
tunnel between the center and the SiSET. The resulting change in /set can be understood 
by considering the donor-SET hybrid system as analogous to a double quantum dot in the 
parallel configuration [15] . We find that the amplitudes in a pattern of Coulomb peaks 
depend on the pulse duration and duty cycle, relative to the emission and capture rates for 
tunneling from or onto the donor. Employing a rate equation model, we are able to extract 
the electron tunneling rate for a specific charge center. We observe fast tunneling rates of 
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FIG. 1: Scanning electron micrograph of a hybrid device. The P donors are implanted close to an 
induced SiSET (grey dashed square). The SiSET is formed by two gate controlled (Vbi,Vb2) tunnel 
junctions and the overlapping top-gate. Panel (b) displays a close up of the stability diagram near 
the charge transition of a center with large charge transfer signal Aq ~ 0.6e. A large gate voltage 
scan is shown in (c), where various charge transitions from multiple centers are visible. 

3000 s -1 for the tunnel event to load the center and 1000 s -1 for the reverse process, despite 
the detection bandwidth for the SET (DC) readout being limited to 200 Hz. 

Figure [I] (a) shows the device fabricated on a high purity intrinsic silicon wafer (> 
10 kficm), with the implantation sites (grey dots) located next to an induced SiSET [14, 16J. 
In the active device region a high-quality, 5 nm thick silicon oxide is grown by dry thermal 
oxidation, yielding a very low density of interface traps ~ 10 10 /eV/cm 2 near the conduction 
band edge p2]. Underneath this oxide ohmic contacts are provided ([P]=5xl0 19 cm -3 ). In a 
first electron beam lithography (EBL) step, with subsequent development and evaporation, 
Ti(15 nm)/Pt(65 nm) alignment markers are formed for a high precision (< 20 nm) realign- 
ment of subsequent layers. A 90 x 90 nm 2 aperture is opened in the PMMA resist, acting as 
mask for the 31 P donors, which are implanted with an acceleration voltage of 14 keV and at 
a fluence of 2.5 x 10 11 cm -2 , resulting in a total of 20 ± 4 31 P donors in this region. After a 
rapid thermal anneal (1000 °C, 5 s) to repair the implantation damage, the Al donor control 
gate as well as the Al barrier gates of the SiSET are patterned. The surface of these gates 
is oxidized by an O2 plasma ash for 4 min at 180°C, resulting in a ~ 5 nm thick Al^O^ 
insulating layer [IB]. An Al top-gate, overlaying the barriers and the source-drain regions, 
is formed in the last EBL step. This process results in a hybrid quantum system with a few 
31 P donors in close vicinity to a SiSET. The sample is operated in a dilution refrigerator 
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at an electron temperature ~ 200 mK. The source, drain, as well as the SiSET control 
gates are connected to the room temperature electronics via Cu powder filters with a cut 
off frequency ~ 1 GHz. The donor control gate is connected via a high-bandwidth line to 
apply high frequency pulses (/mb ~ 500 kHz, limited by resistive voltage dividers at room 
temperature). Its voltage Vb is the sum of a constant component plus a rectangular wave 
for pulsed voltage spectroscopy. The SiSET DC source-drain current is measured using a 
current amplifier with 200 Hz bandwidth and a gain of 10 10 V/A. 

The capacitively- and tunnel-coupled donor and the SiSET island effectively form a double 
quantum dot in parallel configuration [15J. Both series and parallel configurations result in 
a hexagonal stability diagram, but in the series configuration, transport only occurs at the 
triple points[19j. In contrast, in the parallel configuration the transport channel is open 
for any gate voltages for which the electrochemical potential of the SiSET, //set, resides in 
the source-drain bias window. As a result, transport occurs along some of the lines that 
connect the triple points, which we call transport-lines in the following. Fig. [T] (b) shows 
these transport-lines in the vicinity of a charge transition, measured with a source-drain 
bias Vsd = 50 /xV. The relevant gate space is defined by the top-gate (Vr) of the SiSET 
and the donor control gate (Vb). When the energy level of the donor is raised with respect 
to /isETj at the charge transition point it becomes favorable to remove an electron. This 
change in the charge configuration (here labeled as D° — > D + transition) acts back on /zset 
and results in a shift of the Coulomb peak lines. The magnitude of the shift in /isET> relative 
to the Coulomb peak spacing, is quantified by the charge transfer signal Aq w 0.6e. Since 
Aq is much larger than the width of the Coulomb peaks, /set is switched from zero to its 
maximum value by changing the occupancy of the charge center. The data in Fig. [I] (b) 
demonstrate the ability to resolve with essentially 100% contrast the charge state of the 
donor, a critical prerequisite for the spin readout method proposed in Ref. [121 

As shown in Fig[l] (c), the measurement of /set as a function of Vp and Vb yields a set of 
Coulomb peaks appearing as tilted lines (due to the cross-capacitance between control gate 
and SET island) that break at the charge transition points. For Vb > -0.6 V the slope of 
the transport-lines decreases, indicating charge accumulation under the donor control gate. 
In this regime we find several small charge transitions with Aq < 0.1, which we interpret 
as the ionization of shallow charge centers to the Si/Si02 interface. At more negative 
voltages, the pattern clears up, showing well-isolated charge transfers with 0.2 < Aq < 0.6 
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FIG. 2: Charge stability diagram for the pulsed voltage spectroscopy, (a) Sketch of the SiSET 
conductance as function of the top gate voltage Vt and the dc component of the donor control gate 
Vdo- Drawn as blue (green) lines are the positions of the Coulomb peaks when the added pulse 
voltage Vp=Vl=0 V (Vp=Vh). The dashed lines are guides to the eye to indicate a slice of the 
hexagonal charge stability diagram characteristic for a double dot system. When the frequency 
of the square wave /p for switching between Vl and Vh is much smaller than the tunneling rates 
r c and r e , no current is observed between the two superimposed stability diagrams (e.g. panel 
(b), where /p = 61.3 Hz with duty cycle d = 0.5), because the charge configuration can follow the 
equilibrium state. In contrast, when 2irfp 3> r c , T e a non-equilibrium charge configuration can 
be observed resulting in current along the orange lines. Panel (c) displays /set for /p = 5.12 kHz, 
where additional current is visible in this area. Panel (d) shows the schematic for the pulsed voltage 
imposing a charge transition, indicating the various relevant times for the rate equation model as 
described in the text. 

in agreement with the predicted values for electrons tunneling into the SET island from 
donors ~ 30 — 50 nm away [12] , and similar to the values observed in [20] for a charge center 
near Al- and Si-SETs. This part of the stability diagram is stable and reproducible upon 
thermal cycling to room temperature. 
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We stress that the parallel geometrical configuration of our hybrid device impedes di- 
rect spectroscopy of the charge center coupled to the SiSET. It is therefore impossible to 
distinguish a donor from e.g. an interface trap on the basis of its energy levels structure 
[2T] . The number of charge transitions observed for Vd < -0.6 V is compatible with the 
number of donors expected to be found within 30 to 50 nm from the SET island, given the 
P implant fluence. Furthermore, we note that the charge transitions in this regime typically 
group in pairs, agreeing in Aq and slope in the Vr-Vb gate space, again compatible with the 
observations of P donors with two charge transition levels expected [H] 122]- However, the 
unambiguous identification of the charge center remains a quest for spin readout in com- 
bination with magnetic resonance techniques [12] . In this letter, the main focus is on the 
study of tunnel rates between a SiSET and a charge center, whose precise nature does not 
affect the results. 

We measure the tunnel rates by superimposing on the DC voltage of the donor control gate 
(Vbo) a rectangular wave with frequency fp, duty cycle d, and amplitude Vp (cf. Fig. |2](d)). 
If fp is slow compared to the (de) charging rate of the center, we record two stability diagrams 
(blue and green in Fig. [2] (a)), offset by Vp on the horizontal axis when plotted vs. the DC gate 
voltages Vbo an d Vp. These arise because any point on the diagram probes the average /set 
for the combination of the gate voltages (Vbo + V, Vr) and (Vbo + Vp, Vp). Conversely, 
if fp is faster than the electron tunnel rate to/from the charge center, we find /set 7^ 
at gate configurations where transport would be otherwise suppressed, which we call non- 
equilibrium transport-lines in the following (cf. orange lines in Fig. [2]), in addition to the 
pure shifting of the pattern along the Vbo-axis. These lines arise because the charge center 
retains its configuration for the time-span determined by the tunneling time, even while its 
chemical potential crosses the charge transition point. To be specific, at Ib (cf. Fig. [2] (a)) 
no current is expected for the D + configuration, the equilibrium state at Vp = 0. When 
the additional voltage pulse is in the high state, the chemical potential is pushed over the 
charge transition point, into the region where D° is the equilibrium configuration. If an 
electron is captured (and D° is occupied) during this time, immediately after Vp is brought 
back to zero we will find the D° state at a gate configuration where a transport-line is 
present. Thus, /set 7^ when Vp = 0, until the electron tunnels out again. Observing a DC 
current in the orange shaded area around Ib indicates that the system is able to maintain 
a non-equilibrium state for a time comparable to the pulse duration, i.e. the tunnel rate 
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is comparable to the pulsing frequency fp. Similar arguments hold for Iq and Id, whereas 
the bias line around I a is not altered, because the voltage pulse does not cause a charge 
transition of the center. Note that this symmetry break stems from the asymmetric pulsing 
between V and Vp. The low-frequency limit is shown in Fig. [2] (b), where fp = 61.3 Hz, 
and only a horizontally-shifted duplicate of the Coulomb peaks pattern is observed. In 
contrast, the data in Fig. [2] (c) illustrate the high-frequency limit fp = 5.12 kHz, where we 
find /set 7^ at the location of non-equilibrium transport lines. 

Quantitatively, the DC value of /set at the non-equilibrium current peaks can be under- 
stood within a rate equation model. When we pulse the chemical potential of the charge 
center over the charge transition level, the current state D + /° will either persist, because the 
stable state is reached, or change to the opposite state with the corresponding capture (r c ) or 
emission (r e ) rate. The probability to find at the point j in Fig. [2] (d) the D + state occupied 
is Pj(D + ) = Pj(D + )exp(— T c th), because during r H the system tends towards D°. Addi- 
tionally, Pj.(D°) = P T (D°) + P T (D+)(1 - exp(-r c Ttf)), because P;(D°) is increased during 
tjj. The same arguments hold for the inverse direction yielding Pj(D°) = Pj(D°)exp(— T c tl) 
and Pj(D + ) = Pj(D + ) + Pj_(D°)(l — exp(— T e tl)). To determine the four separate time dura- 
tions, corresponding to occupation of either D° or D + , for both values or Vp, which we label 
Tl/h(D° /D + )) we express the probabilities Pjyj(D°/D + ) as a function of T e and T c . We in- 
tegrate these probabilities, including the time evolution of the occupation of the charge state 
over the pulse length tl/h, to obtain the average time, finding D°/D + during Vp = Vq/Vh- 
The result is the four times of interest: T L (D°) = (l/T e )S, T L (D + ) = t l — (l/T e )S, 
T H (D+) = (1/F C )S, and T H (D°) = r H - (1/F C )S, where S = ^^fig^r^?" - 

The times obtained in this manner are proportional to /set- Since every transport-line 
has an individual current amplitude, we analyze the peak ratios i ^ Ib arid jJ° Id which are 
equal to the ratios TH+ 7^(£>o) and T ^^"j^^ D+ ^ , respectively. Although pulsing is performed 
parallel to the control gate axis, it is possible to compare current amplitudes from a cut along 
the top-gate axis, because the current amplitude does not vary significantly along specific 
transport-lines. Figure[3](a) displays the peak ratios for a duty cycle of d = 0.5 as function of 
fp. The ratios t J+ Ib (black squares) and jj+ ln (green squares) are obtained from data like 
Fig. (2](b) or (c) and both show a quantitative agreement with the model using T e = 1000 s _1 
and T c = 3000 s _1 over the entire frequency range fp. Fig. [3] (b) compares the model with 
the experimental data for a fixed fp = 61.3 Hz as a function of the duty cycle d. Again, the 
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FIG. 3: Intensity ratios of the Coulomb peaks for points in the charge stability diagram (cf. 
Fig. [2] (a)) and time ratios originating from the rate equation model (dashed lines). Panel (a) 
shows the peak ratios as functions of the pulse frequency fp for a duty cycle of 0.5. The dashed 
lines represent the model using a capture rate r c = 3000 s _1 and an emission rate T e = 1000 s . 
Panels (b) and (c) show the same ratios as function of the duty cycle for a fixed fp of 61.3 Hz 
and 613 Hz, respectively, where experiment and model is in good agreement. Additionally, the red 
filled circles are the Coulomb peak ratios outside the region of the charge transition. They show 
good agreement with the ratios of tjj and tl depending on the duty cicle d. 

data is described well by the model using the same capture and emission rates. At d « 
and d ~ 1, the peak ratios are more difficult to determine due to the low /set for one of 
the contributions, explaining the deviations from the model. For comparison, the duty cycle 
(red circles in Fig. [3] (b),(c)) is recovered from the spectra independently by analyzing the 
ratio jJ+ If ; showing good agreement with the duty cycle applied. Fig. |3j (c) shows the same 
plot as Fig. [3] (b) for a higher fp = 613 Hz, again in good agreement with the model. 

An estimate of the distance between the charge center and the SET island can be obtained 
from the capacitive modeling of the charge transfer signal Aq, as shown in Ref. Q2J For the 
specific geometry of the device measured here, we find that Aq ~ 0.5e corresponds to a 
distance ~ 40 nm. We use ISE-TCAD to calculate the profile of the conduction band 
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between donor and SET when the D° state is aligned with /xset, and from this, the area of 
the tunnel barrier. A WKB calculation of the tunnel rate yields T ~ 10 4 s _1 , in reasonable 
agreement with experimental findings. 

In summary, we demonstrated and analyzed the tunneling of electrons in a hybrid device 
consisting of 31 P donors implanted next to a gate-induced SiSET. We showed that the 
changes in the surrounding charge configuration can be sensitively detected by the SET, 
and the mutual coupling fulfills the requirements necessary for spin readout as proposed 
in Ref. We further demonstrated a technique to determine the tunnel rate of the 

center investigated, and this technique is applicable even when this tunnel rate exceeds the 
bandwidth of the detection SET. We also provide a quantitative tunnel rate model that 
agrees with the experimental findings. This experimental and theoretical toolbox paves the 
way to the use of spin-dependent electron tunneling as a readout method for single spins in 
silicon. 
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